NASA TECHNICAL NOTE 


CO 

CO 

U-l 

^r 


<c 

on 



NASA LN P-4538 

d. / 



LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
KIRTLAND AFB, N MEX 


BOUNDARY CORRECTIONS FOR A 
THREE-COIL CONDUCTIVITY/ VELOCITY 
PLASMA PROBE 


by Edward W. Vendell 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


APRIL 1968 


tech libraby kafb, nm 



tech library kafb, nm 



0131127 

JNASA TN D-4538 


BOUNDARY CORRECTIONS FOR A THREE-COIL 
CONDUCTIVITY/VELOCITY PLASMA PROBE 


By Edward W. Vendell 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


or sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 3 

A THREE-COIL CONDUCTIVITY/VELOCITY PROBE 4 

FIELD EQUATIONS FOR A CYLINDRICAL CONDUCTING FLUID 10 

Electrical Conductivity 10 

Product of Conductivity and Velocity or aUjj 15 

A METHOD FOR COMPUTING CONDUCTIVITY AND VELOCITY PROFILES OF 

AXISYMMETRIC PLASMA JETS 19 

CONCLUDING REMARKS 1 2 6 

REFERENCES 28 



BOUNDARY CORRECTIONS FOR A THREE-COIL 


CONDUCTIVITY/VELOCITY PLASMA PROBE* 
By Edward W . Vendell 
Ames Research Center 


SUMMARY 


A three-coil plasma probe that measures both the electrical conductivity 
and velocity of laboratory plasmas having low magnetic Reynolds numbers has 
been developed and tested by Rossow and Posch. As a first approximation, it 
was assumed that the plasma boundary was far from and much larger than the 
probe. At the suggestion of V. J. Rossow, the present work was undertaken to 
extend the previous theory by deriving factors that correct for the presence 
of cylindrical boundaries. As a check on this numerical work, several com- 
puted values were compared with electrical conductivity measured in cylinders 
of acid. Since the agreement was satisfactory, the boundary correction fac- 
tors were used to reduce data taken as the probe was swept through an argon 
plasma generated by a constricted-arc wind tunnel . These resultant profiles 
represent local conductivity and velocity values for a free plasma jet having 
a cylindrical boundary. It was found that the raw data underestimates 
conductivity and overestimates velocity. 


INTRODUCTION 


The measurement of the electrical conductivity and velocity of high- 
temperature laboratory plasmas has necessitated the modification of tradi- 
tional immersible transducers. Instruments, such as pitot tubes and 
thermocouples, must be designed to withstand high heat-flux rates and to func- 
tion properly in an ionized environment without greatly perturbing the medium. 
In the past, several external measurement techniques have been developed to 
determine electron and ion temperatures, particle number densities, and the 
total electron collision frequency. These properties have been used to infer 
values of the transport coefficients by means of appropriate formulae. 

Measurements within the plasma stream are desirable for checking the 
external methods and for obtaining, if possible, local values. Consequently, 
many conductivity instruments have been conceived (see, e .g . , refs. l-l6), a 
few of which are immersible (see refs. 4, 5, 11, and 16) . A similar situation 
exists for velocity measurement techniques (see, e.g., refs. 12, l4, and 


*This report was prepared while the author was on leave-of -absence from 
Utah State University and was submitted to Oklahoma State University in par- 
tial fulfillment of the requirements for the degree of Doctor of Philosophy 
in Mechanical Engineering, May 1967 * 



16-20) . For a brief review and comparison of conductivity and velocity 
techniques, the reader is referred to reference 21. 

The present work is based on a design by Rossow and Posch (ref. l6) of 
an immersible three-coil electrical conduct ivity/velocity probe that repre- 
sents a significant improvement upon previous methods because it minimizes 
the heat-flux sensitivity (refs, k- and 5) and large flow-perturbation problems 
(ref. k) of other designs and because data taken with this instrument can be 
reduced to obtain local values of conductivity and velocity. Briefly, this 
device consists of three small coils, one primary and two secondary. The 
alternating current in the primary coil creates an oscillating magnetic 
dipole field. The secondary coils are located in spatial positions such that 
the voltage induced on one coil is linearly related to the conductivity while 
the voltage on the other coil is linearly related to the conductivity- 
velocity product. Thus the data taken as the probe sweeps across a diameter 
of a plasma jet may be used to determine both conductivity and velocity 
profiles . 

In the theoretical analysis of most magnetofluidmechanic problems, the 
magnetic Reynolds number, R m = apUZ, indicates the relative magnitudes between 
the impressed and induced magnetic fields. In this definition a is the 
electrical conductivity, \i is the magnetic permeability, U is the speed, 
and Z is some characteristic length. If R m is small, reference 22 shows 
that the induced magnetic field is also small compared to the impressed field. 
Then the analysis can be greatly simplified by the use of a power series 
expansion in R m since only first-order terms need be considered. This 
approach was used in the analysis of Rossow and Posch because R m is small 
for most plasmas generated by electric arcs; the probe design of reference 1 6 
was tested in a constricted-arc wind tunnel where R m ranged from 10“ 4 to 
10“ Therefore, the neglect of terms of order R m 2 was reasonable and 
these terms will also be omitted in the present analysis. 

The theory of Rossow and Posch also assumed an unbounded plasma having 
uniform electrical conductivity and velocity. The principal purpose of this 
paper is to present theoretical modifications which will remove these restric- 
tions for the case of a free plasma jet having a cylindrical boundary. 

Accordingly, modifications of the unbounded field theory that are neces- 
sary whenever the three-coil probe nears a plasma boundary are derived and 
computed numerically. Initially, the conductivity a and velocity U are 
assumed to be constant inside the jet of radius R and zero elsewhere. From 
this idealized model, the modifications are derived and presented graphically 
in the form of boundary correction factors. These theoretical correction 
factors are compared with experimental results obtained by immersing a three- 
coil probe in a deep plastic cylinder filled with an electrolyte. 

Then the correction factors are used to develop a method for computing 
the local conductivity and velocity of a plasma stream from data which varied 
continuously across a low-density argon plasma jet generated by a constricted- 
arc wind tunnel. The method is then applied to a typical data record and the 
resulting conductivity and velocity profiles are presented and discussed. 
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SYMBOLS 


A 

B 

B 

~P 

St'SlI 

C S 

c Y 

D 

E 

e 

H 

I P 

J 

n p 

n E 

n Y 

P 

P 

R 


r S 

r Y 

S 


vector potential function 
magnetic induction vector 
primary dipole magnetic induction vector 
first-order magnetic induction vectors 
conductivity calibration constant 
velocity calibration constant 
electric displacement vector 
electric field intensity vector 
electron charge 

magnetic field intensity vector 
total current, primary coil 
current density vector 
first-order current density vectors 
number of turns, primary coil 
number of turns, E-coil 
number of turns, Y-coil 

subscript referring to primary-coil parameters 

dimensionless radial distance to location of primary coil (P = p/R) 

radius of a cylindrical region 

magnetic Reynolds number 

distance in spherical coordinates 

characteristic radius, primary coil 

characteristic radius. E-coil 

characteristic radius, Y-coil 

dimensionless coil spacing parameter (S = y^/R) 
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subscript refers to first-order fields caused by a 

velocity vector 

axial velocity component 

dimensionless Cartesian coordinates 

Cartesian coordinates 

radial distance from primary coil to either secondary coil 

axial coordinate 

dielectric constant 

boundary correction factor for a 

boundary correction factor for crtJ 

azimuthal coordinate 

magnetic permeability 

radial distance in cylindrical coordinates 
conductivity- coil parameters 
electrical conductivity 
conductivity velocity coil parameters 
signal induced on E-coil 

signal induced on E-coil in an infinite medium 
signal induced on Y-coil 

signal induced on Y-coil in an infinite medium 
stream function for current density vector 
frequency of impressed power 
first-order fields caused by crU|| 

A THREE-COIL CONDUCTIVITY/VELOCITY PROBE 


Rossow and Posch (ref. 1 6 ) have developed and tested an immersible, 
three-coil plasma probe that represents a significant improvement upon 
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previous techniques because it is capable of determining local electrical 
conductivity and velocity values without greatly perturbing the flow and 
because it can be used in such a manner that heat-flux sensitivity is neglig- 
ible . Since the present work is based on the probe design of Rossow and 
Posch, selected material from reference 1 6 will be reviewed in this section. 


Y 



Figure 1.- Perturbation of Bp caused by a 
conducting fluid at rest. 


Y 



Figure 2.- Perturbation of Bp caused by a moving 
conducting fluid. 


The basic concepts of the instru- 
ment are illustrated in figures 1 and 
2 which were prepared by superposition 
of the probe onto figures 2 and 3 of 
reference 1 6 . Figure 1 represents the 
principle by which conductivity is 
measured. Three small coils labeled 
P for primary, E for conductivity, 
and Y for velocity, lie in the same 
plane (x = 0) and are mounted on sup- 
port rods in such a way that the axes 
of the E- and P-coils are parallel to 
the z axis while the axis of the 
Y-coil is parallel to the y axis . 

In the absence of a conducting or 
a dielectric medium, the pulsating 
current in the primary coil produces 
magnetic flux which may be approxi- 
mated by an oscillating ideal magnetic 
dipole field. The solid lines labeled 
Bp indicate the shape of the lines of 
force for the primary dipole field at 
an instant in time. Under these cir- 
cumstances, the primary field, having 
no axial component, produces no flux 
linkage at either the E-coil or the 
Y-coil. However, in the presence of 
a conducting, quiescent medium, the 
primary magnetic field induces cur- 
rents which, in turn, create a per- 
turbation magnetic field designated 
b^_. The dotted lines indicate an 

instantaneous position of the lines of 
force for the resultant magnetic field 
which produces flux linkage at the 
E-coil; however, the Y-coil remains 
free of flux linkage because of its 
orientation. 

Similarly, in figure 2, the 
dotted lines indicate a particular 
instantaneous position of the lines 
of force for the resultant magnetic 
field which is the sum of the primary 
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field Bp and a perturbation field b|| caused by the motion of a conducting 
fluid across the primary lines of force. As the sketch shows , the resultant 
field produces flux linkage at the Y-coil while the E-coil flux linkage is 
zero for such a disturbance. 

Since the perturbation fluxes at the two secondary coils are time depen- 
dent, Faraday's law of induction implies that the potential output of the 
E-coil will be proportional to the magnitude of the z component of the field 
perturbation caused by the conductivity a, and the potential output of the 
Y-coil will be proportional to the magnitude of the y component of the field 
perturbation caused by the product of the conductivity a and velocity Un 
or aU||. 

The theoretical analysis of Rossow and Posch begins with Maxwell's 
equations 


the simplified Ohm's law 


V • B =• 0 , V • E = 0 

V X B = [iJ 

SB 

V X E = - — 


J = ct(E + U X B) 

the conservation of charge equation for a neutral plasma 


V • J = 0 

and the Coulomb or transverse-gage condition 

V • A = 0 


( 1 ) 

( 2 ) 

(3) 

W 

(5) 

( 6 ) 


In these equations E is the electric field intensity vector, B is the mag- 
netic induction vector, J is the current density vector, U is the plasma 
velocity vector, a is the electrical conductivity, \i is the permeability, 
and A is a vector potential function such that 

B = v x A (7) 

To obtain a solution for these equations, Rossow and Posch imposed the 
following restrictions; 


1. The flow field is unbounded. 

2. The electrical conductivity and velocity are taken to be constant 
over the entire flow field. 
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3. The only applied field is B p which is the magnetic field created 
hy the oscillating current in the primary coil. 

4. Since the magnitude of Bp is less than 10 gauss , the assumption 
of scalar electrical conductivity is justifiable . (in the presence of large 
applied fields, the conductivity of a plasma assumes a tensor ial form.) 

5. The small magnitude of Bp also justifies the representation of the 
resultant magnetic field, B, as a power series expansion in the magnetic 
Reynolds number, R m . As mentioned previously, typical values of R m for low 
density constricted-arc wind tunnels range from 10~ 4 to 10~ 1 so that neglect- 
ing terms of order R m 2 is reasonable. (See ref. 22 for details on series 
expansions in the magnetic Reynolds number.) 

6. Perturbations of the stream caused by the presence of the probe are 
neglected. 

7- In equation (2) the displacement current term is neglected. 

Omission of this term is justified in reference 21. 

As a zeroth-order approximation for B, Rossow and Posch used the primary 
dipole field B that would be produced by an idealized primary coil in 

free space. The vector potential A^ for such a field is well known and 
appears on page 237 of reference 23 as 

A p (x,y,z,t) (8) 

where m is the magnetic dipole moment of the idealized primary coil located 
at (x,y,z) = (0,0,0), and r = (x 2 + y 2 + z 2 ) 1 / 2 . Thus, equation ( 7 ) yields 
the solution 

~P = ~ X ~P ( 9 ) 

Since Bp is time dependent, equation (3) requires the existence of an 
associated electric field, denoted by subscript t, such that 

dB D 

v x E t = - ^ (10) 

and this equation may be solved for . 

Equation (^) suggests that the total current may be considered the sum 
of two components, the first of which is 

J t = crE t (ll) 


and the second is 
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j|| = ^ X Bp = trU||k X Bp 


( 12 ) 


where U = U||k, as in figure 2 with the unit vectors i ; k directed in 

the positive x, y, z directions, respectively. The first component, J is 
caused by the application of an electric field, to a stationary conducting 
fluid. The second component, J|| , arises from the motion U||k of a conducting 
fluid across lines of force, B^. 

Rossow and Posch continued the analysis by using equation (2) to solve 
for the first-order perturbation fields b^. and b|| that accompany and 

J||. Thus, the resultant magnetic field, including first-order terms only, is 

5 = 5p + St + 5|| ( 13 ) 


This equation, together with equations (2), (3), and (*0 , could be used to 
solve for currents and perturbation fields of higher order in R^. However , 
as previously mentioned, the terms of higher order may be neglected in the 
present case because is small. 


For future reference, the following solutions are listed from 
reference 1 6 : 


B^ = -m cos co t 


. 3 xz . 3 J z fe / 3 z 2 n 
h — + J — + ”3 l — " 1 
r ~ r \ r 


Et = -com sin cotV X ( r* 


J-j- =■ -a com sin cot 


/V\ ^ . /x 
iW + £1^3 


J,, = 3U||m cos cot 


. yz . xz 

i 5 " £ 5 

. r r . 


t>4- = 


a^com 


sin cot 


. xz . yz 

i s 


r 2 + z 2 


an u n m 

b 1 1 = — - — cos cot 


1 - 


3£ 
r 2 J 


. x , . y , z 

1 3^"J 3 + k q 

r 3 ~ r 3 r . 


m = (n p r p 2 )M.Ip/4 


(d) 

(15) 

( 16 ) 

(17) 

( 18 ) 

(19) 

( 20 ) 
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where I p is the peak current supplied at frequency u> to the n p -turn 
primary coil whose characteristic radius is r . 

Equations (l4) through (20), plus Faraday’s law of induction, suggest 
that locations and orientations of the E- and Y-coils may he selected so as 
to eliminate unwanted magnetic flux linkage. Faraday’s law states that the 
potential induced on a secondary anywhere in the flow field may he 

approximated hy 


0 S 


“TtnoT 


s- 1 - s 


lr « 


Ss) 


( 21 ) 


where r s is the characteristic radius of the n s -turn secondary, B is the 
resultant field evaluated at the center of the coil, and N s I s a unit vector 
parallel to the axis of the coil. Thus, if the center of the E-coil is 
located at (x,y,z) = (0 ,y^ ,y^,/\/2) and is oriented so that Ng = k, equa- 
tions (l8) and (21) may he combined to obtain 


2 Jl - L p u 

$ Eoo = ( n l5£ 2 ) ( n p r p )(l^ w ) 2 GOS 

6 s/ 3/2 y 2 

where 00 refers to an unbounded medium and t is time. Similarly, 
center of the Y-coil is located at (x,y,z) = (0,-y^,0) so that Ny 

the potential induced on the Y-coil in an unbounded medium will be 


( 22 ) 

if the 
j, then 


p^unlpcrUj, 

<t> Yoo = - ( n y r Y 2 ) ( n P r p 2 ) i sin wt ( 2 3) 

8y£ 

Thus, equation (22) indicates that is linearly related to a and 

equation ( 23 ) displays the linear relationship between <4^ and crU||. 


Either one of the two calibration methods given in reference 1 6 may be 
used to present equations (22) and ( 23 ) in the equivalent form 


•&» 


a 


(24) 


and 


^00 

aU tl 


(25) 


where the calibration constants Cy> and Cy are functions of the various 

probe parameters as well as the associated electronic circuitry. It is con- 
venient to compute and Cy so that and represent peak-to-peak 

values as determined from an oscilloscope data record. 


In the next section equations (2^) and (25) will be further modified by 
the introduction of factors that correct for the presence of a cylindrical 
boundary. 
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FIELD EQUATIONS FOR A CYLINDRICAL CONDUCTING FLUID 


As mentioned above, the theory of Rossow and Posch (l6) decoupled the 
perturbation magnetic fields by neglecting terms of order R^ 2 , assumed a 
and Uj| to be constant and uniform, and considered only an unbounded medium. 
This section presents modifications of the previous theory by recognizing the 
presence of a cylindrical boundary of radius R such that the plasma occupies 
the region 

x 2 + y 2 < R 2 ; -oo < z < oo 

Also, the position of the probe will not be restricted to the center of the 
cylindrical region. Furthermore, for this part of the analysis it is assumed 
that a and Ujj are constant and uniform in the region occupied by the plasma 
and that these quantities vanish elsewhere. The present analysis will also 
ignore terms of order R m 2 , the displacement current, and relativistic 
effects. The only applied field is Bp, which is caused by the oscillating 

current in the primary coil, and this field, being of the order 10~ 2 gauss, 
is so small that the assumption of scalar electrical conductivity is justi- 
fied. Finally, as above, perturbations of the stream caused by the presence 
of the probe are neglected. 


Electrical Conductivity 

The analysis begins with a consideration of the basic equations which 
are identical, except for the boundary conditions, to equations (l) through 

( 7 ). 


Because 0 and U|j 
ditions on the current 

vanish 

density 

outside the cylinder, the two boundary 
vector are 

con- 


11 

10 

whenever y? + y 2 > R 2 

(26) 

J 

11 

0 

for x 2 + y 2 = R 2 

(27) 


where the unit vector, N_l, is perpendicular to the cylindrical boundary. 

If equation (7) is substituted into equation (2), the result will be 

V x V x A = pJ ( 28 ) 

and if this expression is expanded so that equation (6) may be applied, the 
final result will be 

v 2 a = -\±J (29) 


10 



where it is understood that the Laplacian operates on each rectangular com- 
ponent of A. Equation (29) represents a condensation of the three Maxwell 
equations. Therefore , the problem reduces to the solution of equations ( 4 ) 
and (29) subject to the boundary conditions given by equations (26) and (27) . 

Equation (29) has a unique solution and it may be solved by finding a 
Green* s function or by an equivalent technique, the method of images (see, 
e.g., ch. 2 , ref. 2 * 0 . The latter method was chosen because it appears to be 
a simpler approach and it makes use of equations (l 4 ) through ( 20 ) . 

It is also convenient to use an analogy based on the steady, two- 
dimensional flow of an incompressible fluid for ^which there exists a potential 
function ¥, called the stream function, such that the velocity field U is 
given by 


U = V X k! 


(30) 


An important property of the stream function is the fact that lines of 
constant 1 are streamlines. 

This analogy may be used to represent any two-dimensional current density 
vector field as 


J = v x k\|/ 


(31) 


with the convenient property that current paths or loops coincide with lines 
of constant >|r. 


Returning to the method of images, con- 
sider the geometry of figure 3, where the medium 
•is assumed to be unbounded, at rest, and of uni- 
form and constant a. If the real coil is 
located at (x,y,z) = (a,0,0), an image system 
of strength A at (x,y,z) - (-a, 0,0) must be 
found such that one of the induced current loops 
will coincide with the dashed cylinder and 
thereby satisfy the boundary condition given by 
equation (27). 

Figure 3.- coil and image system The current, J^ r , induced by the real coil 

is given by equation (l 6 ) which, for the 
present case, becomes 



where 

G-fc = crwm sin urt 

r 2 = (x - a ) 2 + y 2 + z 2 
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and 


itr 


m is defined by equation (20 ) . 
is 

*tr = 


Therefore, a real stream function for 
a) 2 + y 2 + z 2 ]" 1/2 (33) 


By trial and error, a stream function 


tti =AG t [(x + a) 2 + y 2 +A 2 z 2 ]” 1/2 (3^) 

vas found for the image system so that the resultant stream function is 

i|r t = -G t {[(x - a) 2 + y 2 + z 2 ]~ 1/2 - A[(x + a) 2 + y 2 + A 2 z 2 ]' 1/2 } (35) 

and the total induced current becomes J^_ = v X This expression for ^ 

is particularly useful because the surface defined by the condition ^ = 0 
happens to be the cylindrical surface 


x - a 


A + 1 
.A 2 _ i/J 


v2 


+ y 2 » 


2a A 




(36) 


so that the current loops corresponding to ^ = 0 meet the requirement 
specified by equation ( 27 ). 


y 



Figure 4 illustrates the locations of the 
real primary and secondary E-coils after the y 
axis has been translated a distance 
L = a(A 2 + l)/(A 2 - l) , that is, to the center 
of the jet . Since L - a = p and 
R = 2aA / (A 2 - l) , a and A may be solved as 
functions of R and p with the result that 


Figure 4.- Position of coils inside 
a cylindrical region. 

a = (R 2 - P 2 )/2p 

(37) 


A = R/p 

(38) 


If these expressions are substituted into equation (35) and if the y axis 
is translated a distance L, it follows that 


*t 


= -G-fc[( x + p 2 + y 2 + z 2 ) ^ 2 - R(px + R 2 + p 2 y 2 + R 2 z 2 ) -1//2 ] 


(39) 


Therefore , the final expression for the resultant current density vector is 
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I 



I 


Jt = v X k\|r t 


- -GjA 1 


Rp 2 y 


■(px -f-R 2 + p - 2 y 2 + R 2 z 2 ) 3 / 2 (x+ p 2 + y 2 + Z 2 ) 3 / 2 J 


Rp(px + R 2 ) 


x + p 


(px + R 22 h- p^y 2 + R 2 z 2 ) 3//2 (x+ p 2 + y 2 + z 2 ) 3/2 J 


(4o) 


Equation (2) may now be used to obtain b^ from the relation 

V X b t = p.J t (4l) 

by use of equation (ll), of reference 23- That is, the perturbation field, 
b^ , at any field point (x*,y f ,z T ) in the cylindrical plasma is given as 




n/r 2 -- 2 


5t(x' , y ' , z ' ,t) = 


■y 


bit 


J t (x,y,z,t) X Vr -1 dx dy dz 


-R 




(te) 


where r = (x - x f2 + y - y 12 + z - z ,2 ) 1//2 . Note that the current density 
vector J, does not vanish outside the cylindrical region as required by 
equation (26). However, this requirement is satisfied by equation (42) 
because the limits of integration do not extend beyond the cylindrical bound- 
ary so that the mathematical current loops outside the cylinder cannot con- 
tribute to b^. Since the axis of the Z-coil is parallel to k, only the 

z component of b^ will be used in Faraday r s law, equation (2l). Since 
(x T ,y T ,z ! ) = (-p,y^,y^/ N /^) , the z component of b^ must be 


■^oo pE p^R 2 -y 2 


J zt 




J-oo^-R 


Rf 


py(y - y 2 ) + (px + r 2 ) (x+ p) 

7; — 2.3/ 2 


X + P y- y s 2 + Z- -j=J (px+ R 2 W p 2 y 2 + R 2 z 2 


3/2 


y(y- y s ) + (x+ p) ; 


- 2 V 3/ 2 


|_(^x+ p +y-y s +z- — 


■) (x + p 2 + y 2 + 



dx dy dz 


(43) 
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If the variables x, y, z and parameters p, y^, in the last equation are 
made dimensionless with respect to the jet radius R, 





X+P +Y- 


P 2 Y(Y - S) + P(X + P) (PX + 1) 
S' 2 + Z -D 2 j 3/2 (PX+ 1 2 + P 2 Y 2 + 



Y(Y- S) + (X+ P) ; 




+ P +Y - S + Z - D 


rc 


2 2 2 

X + P 2 + Y + Z J J 


ax <3Y dZ 


m 


where X - x/R, Y - y/R? Z - z/R, P - p/R, S - y^/R, and D = S/\[2. Substitu- 
tion of this expression into Faraday's law yields 




-(n 2 r 2 2 )(n p r p 2 )(pu))‘ 


la cos wt 

* -lg B ’ 





dX dY dZ 


(^5) 


where F^ is used to denote the complicated integrand of the preceding 
expression. In order to eliminate the probe and coil characteristics, 
equation (45) is made dimensionless through division by equation (22) to 
obtain the desired correction factor 


= e„(s,p) = - 




l6rt 


-1 


on/i-Y £ 


F s dX dY aZ 


(46) 


For a given probe, the conductivity correction factor represents 

the ratio of the E-coil signal that would be induced in a finite cylindrical 
region having a constant conductivity to the signal that would be induced in 
an infinite region having the same constant conductivity. 

As might be anticipated, the integrals of equation (46) cannot be evalu- 
ated readily in closed form. Although one integration could be accomplished, 
it was more convenient to evaluate the integrals numerically rather than deal 
with the elliptic integrals resulting from the analytical integration. There- 
fore, equation (46) was computed using a numerical 10-point Gaussian- 
quadrature program written in Fortran IV and executed on an IBM 7090/7094 
digital computer. (See ref. 21 for further details.) The results of the com- 
putation are presented in figure 5 which is a plot of ©£ vs . P = p/R for 
four practical values of the coil spacing parameter S = y^/ R; the dashed por- 
tions represent reasonable extrapolations of the computed curves. Additional 
points were not computed because each calculation of ©£ consumed about one- 
half hour of computer time. However, the data reduction method presented in 
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1.0 



Figure 5 -~ Conductivity correction factor. 


a subsequent section shows that satis- 
factory results can be achieved with 
the curves of figure 5- Values, 
rounded to two significant figures, are 
listed for convenience in table I. 

The basic conductivity formula 
(eq. (24)) may now be modified to 
account for the presence of a cylindri- 
cal boundary by introducing the correc- 
tion factor with the result that 


a - 

0e(S,P) 

where C^ is a calibration constant to 
be determined by the methods outlined 
in reference l6 . The last expression 
is still limited to the case of con- 
stant electrical conductivity. 



TABLE I.- CONDUCTIVITY CORRECTION FACTORS, 


P = p/R 





0.05 

1 0 
i'o 

0 

OJ 

0 

0.00 

0.95 

0.89 

0.78 

.4o 

.94 

.87 

.74 

.60 

•92 

.83 

.66 

.80 

.85 

.69 

.44 

.90 

•71 

.47 

.18 


S = y^/R 

0.30 

0.80 

0.93 

1.00 

0 .68 
.62 
•51 
.25 
.044 

0.l4 

0.012 

- 0.05 

— 

— 
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Product of Conductivity and Velocity or o\J\\ 

Another application of the method of images yields a aU^ stream 
function \|r M such that 


'I',, 1 1 z I 


( X+ P 


.-1-3/2 r. 

J " R L 


+ y + z‘ 


- R^l px+ R^ + p 2 y 2 + R^z' 




r /2 } 


(w) 


where G|| = aU|jm cos u)t, and the y axis has been translated to the center 
of the stream as in figure 4. Using the stream function analogy explained 
above, one can derive the following expression for the current density vector: 
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= -G ! ( z||^3R 3 pi s y ( px+R^ 2 + p 2 y 2 + P 2 z 2 ) " 5/2 - 3y(x+p“ + y 2 + z 2 ) ” 5/2 ]i 
-[ 3 R 3 P px+R^ 2 ( px+R 22 + p 2 y 2 + R 2 z 2 ) _5/2 - 3(x + p) (x+p 2 + y 2 + z 2 ) “ 5/2 ]j} 


(*9) 


As was the ease with b^, b|| is related to by the triple integral 

fe||(xSySz',t) = J J J J||(x,y,z,t) X Vr -1 dx dy dz ( 50 ) 

where r = (x - x ' 2 + y - y ' 2 + z - z' z ) 1/z . Because the axis of the Y-coil 
is parallel to the y axis, only the j component of b 1 1 will contribute 

to the induced electromotive force. That component, evaluated at 
(x*,y T , 2 i f ) = (-p,-y£> 0 ) and made dimensionless with respect to R, is 


b yn 


-3^G„ 


1 


YZ^ 


^jtR 2 J — OO <J — 1 V 




(X + P + Y + S + Z 2 ) 


3/2 


X 


P 2 

(PX + l 2 + P 2 ^ + Z 2 ) 5 / 2 


1 

(x + P 2 + Y 2 + z 2 ) 5/2 


dX dY dZ 


-3pGn 

4 jtR 2 



( 51 ) 


where Fy represents the integrand. 


Therefore, Faraday's law predicts that the potential 0 y induced at the 
Y-coil will be 


$y = 



3wlph 2 

16 r 2 


sin 


at J J J Fy dX dY dZ 


(52) 


and, if this expression is divided by equation (23) to eliminate probe and 
coil characteristics, 


l 6 





(53) 


i>Y qq 2 


“I 




F Y dX dY dZ 
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The correction factor ©y defined by 
equation (53) is similar to the con- 
ductivity correction factor ©£ 
because, for a given probe, ©y rep- 
resents the ratio of the Y-coil sig- 
nal that would be induced in a finite 
cylindrical region having uniform 
values of a and U| j to the signal 
that would be induced in an infinite 
region having the same uniform values 
of o and U||. For convenience, this 
parameter will hereafter be referred 
to as the velocity correction factor 
even though it applies to the product 
aU, i . 


Figure 6.- Correction factor for aU M 


As with ©g, the integrations 
required by equation (53) were per- 
formed numerically on a digital com- 
puter using a program similar to that 
used for integrations in equation (46). 
The results are presented in figure 6 

and are also listed, to two significant figures, in table II. The curves of 
figure 6 indicate that if S < 0.10, then the deviation of ©y from unity is 
negligible except at the very edge of the cylindrical stream. As anticipated, 
the conductivity correction factors are more important than the velocity 
correction factors. 


TABLE II.- VELOCITY CORRECTION FACTORS, ©y 


P = p/R 

S = y s /R 

0.05 

0.10 

0.20 

0.30 

0.80 

1.00 

0.00 
.60 
.70 
.80 
• 90 
.95 

1.00 

• 99 

• 93 

0.99 
• 9^ 
.75 

0.99 

.98 

.9^ 

.78 

.44 

0.99 

• 97 

• 93 
.84 

.55 

•29 

0.88 

0.73 
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Equation (25) may now be modified to include the velocity correction 
factor: 


aU„ = 




( 5*0 


where Cy is a calibration constant. Equation (5*0 is still limited to a 
conducting fluid having a cylindrical boundary, uniform cr, and uniform U ( j. 

The conductivity correction factors presented in figure 5 were verified 
experimentally using a three-coil probe. (This experiment was performed with 
the aid of Mr. R. E. Posch.) Each of two plexiglas cylinders having differ- 
ent diameters was filled with an electrolyte having a conductivity of 
7^.8 mho/meter. Then the radial variation of the E-coil output was recorded 
for each cylinder. The results were normalized with respect to an experi- 
mentally determined value of (see ref. l6 for details) to obtain 0£. 

are compared with theoretical 
values computed from the curves of 
figure 5- The differences are less 
than 10 percent except near the edge 
where the output of the E-coil was 
very small and was extremely sensi- 
tive to the alinement of the probe* s 
axis parallel to the z axis of the 
cylinder . 

The probe was also placed imme- 
diately adjacent to and outside each 
cylinder wall (in room air) to record 
the magnitude of the signal produced 
by the currents induced in the elec- 
trolyte by the primary dipole field. 
Bp. As indicated by the data points 
at P - 1.025 and P = 1.10, the result 
was that the signals were too small 
to be observed on the oscilloscope at 
its maximum sensitivity setting; that 
is, the induced signals were at least 
ten times smaller than any signals 
obtained inside the cylinders. Now, 
in a free plasma jet it is reasonable to assume that the conductivity 
increases from zero at the edge to a maximum that occurs at or near the 
center line. The result of the experiment indicates that when the probe is 
swept through the jet, the signal induced on the E-coil at a given radial 
position, say p = p , is not signif icantly influenced by the increased con- 
ductivity in the region 0 < p < p Q ; consequently, is primarily a function 

of the variable conductivity in the region p Q < p < R. This observation is 
of fundamental importance to the data reduction technique presented in the 
next section. 


In figure 7 the experimental values of 



Figure 7.- Comparison of theoretical curves with 
experimental measurements in cylinders of 
acid. 
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A METHOD FOR COMPUTING CONDUCTIVITY AND VELOCITY PROFILES 
OF AXISYMMETRIC PLASMA JETS 


In this section a method for computing conductivity and velocity profiles 
is presented and applied to data taken by a three-coil probe in a constricted- 
arc wind tunnel. The method is based on the premise that the continuously 
varying profile can be approximated by a number of steps as suggested in 
figure 8 . It is then assumed that each level can be treated as a cylinder of 
constant conductivity (or constant 



Radius 

Figure 8.- Measured profile subdivided to apply 
boundary corrections . 


aU||J by the theory developed above. 

As noted at the end of the last sec- 
tion , the magnitudes of $£ and $y 
are negligible when the probe is out- 
side a cylinder of conducting fluid. 
Hence , if the analysis of a given pro- 
file is started at the outside bound- 
ary of a jet, the calculations can 
proceed to the center in an explicit 
fashion so that all parameters are 
known as each step is made inward. As 
illustrated in figure 8 , the cylindri- 
cal plasma jet is, as a first approxi- 
mation, subdivided into a finite 
number, n, of concentric cylindrical 
regions each having a different, but 
constant conductivity. The upper part 
of the sketch shows a cross-section of 
the jet and the lower part displays a 
typical ^-trace together with the 
n-step approximation of the actual 
continuously varying conductivity 
profile . 

When the probe is located at the 
radial position p - p x shown in 
figure 8 , the two data points on the 
abscissa axis of figure 7 indicate 
that the influence of the increased 


conductivity in the inner regions has 

a negligible effect on the value of ± . Therefore, 0 ^ ± may be inter- 

preted as the signal which results from placing the probe at a radius p x in 
a cylindrical region of radius p Q containing a fluid of conductivity o 
The signal $2,2 results from placing the probe simultaneously at a radius 
p = p 2 in two cylindrical regions of radii p Q and p ± containing, respec- 
tively, fluids of conductivity cfi and cr 2 - ffi- Similarly, $2,3 cause( T 

by the simultaneous immersion of the probe at a radius p = p 3 in three 
cylindrical regions of radii p Q , p , and p 2 containing, respectively, 
fluids of conductivity a 1, cr 2 - a 1, and cr 3 - cr 2 . The extension of this 
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reasoning to is straightforward and the results may be put into mathe- 

matical form by the following set of n linear equations in n unknowns; 


*2,1 ' ’ PS 


»S, 2 - c| - 5 |) - ® z ( p 7 - ft 




cr, p 3 > 

= cf ®s( > p 7 ) + 
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3 \ . qs q g £3 
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1=0 


C E 


0 ^ 




tl"pT ’ °) ’ n - 


j 


( 55 ) 


In these equations represents the conductivity calibration constant, o Q 

is identically zero, and the argument of is given, respectively, in 

terms of the parameters S = y^/R, and P = p/R. After inserting values for 
©2 from figure these simultaneous equations may be solved for o^> 
l-l, • • • ,? u • 

The same technique may be applied to obtain a aU|| profile and the 
equations that must be solved for an n-step approximation are: 
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where Cy is the velocity calibration constant, [aU||] 0 is identically zero, 
and the values for ©y as functions of S = y^/R and P = p/R are given 
in figure 6. The solutions of equations ( 55 ) and (56) may then he combined 
to obtain an n-step approximation of the U|| profile. As the number of 
subdivisions is increased, subject to the limitations to be discussed below, 
the discontinuous step profiles should approach the actual continuously 
varying profiles. 

The data reduction technique may be divided into three parts as follows: 

1. The edge of the jet may be determined from the E-coil output pro- 
vided that the response time of the probe is known. The response time of a 
probe system may be determined from a consideration of superimposed data 
taken during a two-way sweep through the plasma, stream. Using this technique, 
Rossow and Posch reported a response time of 0 .03 second. Figure 9 i s a 
typical data record obtained by Rossow and Posch with a three-coil probe 

(S = y s /R = 0.234) in a constricted-arc wind tunnel having an arc current 
I a rc = 200 A and the conductivity data (lower trace) indicate that the 
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extreme radius of the jet was about 
9.5 cm. To assess the sensitivity of 
the data reduction method to the mag- 
nitude of the jet radius, profiles 
were computed for three radii: 9.0 cm, 

9.5 cm, and 10.0 cm; the corrected 
center line values of o and U|| for 
the three cases did not differ by more 
than 10 percent. Therefore, by con- 
trast with the measurement technique 
of Donskoi, et al. (ref. 7 ), the 
accuracy of this method is not highly 
dependent on a precise determination of 
the jet radius. 

2. The null signal should be added to the ^ trace and subtracted 

from the trace; reasons for this procedure are discussed below. Then, 

peak-to-peak values of <p£ and <£y may be plotted vs. radial position of the 
instrument; the radial position should be adjusted to compensate for the 
response time of the system. A mean value for the 0y curve is used, thereby 
eliminating the local signal excursions caused by random stream noise. 

3. The next step is to subdivide the cylindrical jet into subregions as 
suggested by figure 8 and to apply equations (55) and (5 6) to the ^ and 0y 
data. It is suggested that this step be repeated several times in order to 
test its accuracy; that is, as a first approximation, use three subregions to 
obtain three-step a and Ujj profiles. Then, as a second approximation, 'use 
four subregions for the computation. Additional approximations obtained by 
increasing the number of subdivisions can be carried out to increase the 
accuracy and definition of the curves. A limitation on the imximum number of 
subdivisions is discussed below. 

The velocity trace of figure 9 indicates a minimum signal at about 
2.8 cm from the center of the oscilloscope screen. This phenomenon is due to 
b || and Bp being l 80 ° out of phase as indicated by equations (ik) and ( 19 ). 
Outside the stream $y nonzero because the actual Bp differs slightly 
from the theoretical Bp- As the probe nears the edge of the stream, the 
magnitude of b^ increases and causes 0 y, which is proportional to 
Bp + b||, to decrease because of the l 8 o°-phase difference. As the probe sweep 
continues into the jet, the magnitude of b|| dominates and $y begins to 
increase. Therefore, when the data are reduced, the $y ordinates should be 
increased by the value of the null signal while the $£ ordinates should be 
decreased by the value of the null signal. 

If the above data reduction technique is applied to the $£ and $y 
signals of figure 9 , then the result is given by figure 10. The dashed 
curves in that figure represent profiles which were not corrected for the 
presence of a cylindrical boundary and are based on the application of 
equations of the form 
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Scope sensitivity 
settings-. 

T ( .02 mv/cm ) 
2 (.05 mv/cm ) 


— Sweep direction 
Figure 9*- Data record. 
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— = Cy and fHli = C Y 
*E3 $ Y 

to center-line values of and $y. The uncorrected U|| profile was 
obtained by dividing the dashed ordinates of the crU^ profile by the dashed 
ordinates of the 0 profile. As anticipated , the conductivity correction 
factor ©jn had a greater influence on the reduction of data than did the 
velocity correction factor ©y. In fact, the example presented in figure 10 




Figure 10.- Electrical conductivity and velocity profiles, "based on the data of figure 9, in a 

free argon plasma jet . 


indicates that the corrections of the aU|| profile are negligible. However, 
because of the conductivity corrections, the corrected profile may 

differ considerably from the uncorrected profile. 

As mentioned above, the choice of subregions is not completely arbitrary 
and it was found that the values of y^ and R impose two important restric- 
tions. The first restriction may be explained by the curves of figure 11. 
Consider a fluid having a cylindrical boundary and, for example, a constant 
conductivity of 888 mho/m. Suppose that an oscilloscope recording is made of 
$£ vs. radial position, p, for this fluid. Then, if the correct values 

are substituted into the formula 

_o_ = °S 

0g(y s /R, p /R) 

for different values of p, the horizontal curve marked Otrue will be the 
result. However, the smallest screen division for many oscilloscopes is 2 mm 
and, with the presence of a slight amount of noise, it is possible to err by 
as much as 1 mm when reading the oscilloscope data record. Therefore, a 
typical scope sensitivity setting of 0.05 mV/cm, a calibration constant of 
3000 mho/m-mV from reference 16, and the above equation were used to compute 
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Figure 11.- Uncertainty resulting from reading 
error . 



Figure 12.- Conductivity profile in an argon 
plasma jet . 


the error curves of figure 11 for two 
different values of the parameter 
y /R. The figure shows that a reading 

error of 1 mm can cause a large error 
in a when y /R = 0.3 and p/R > O.80. 

2j — 

However,, the error is more moderate 
when y^/R = 0.10. 

Thus, if y^/R is small, the 
stream can be subdivided into several 
subregions without risking the intro- 
duction of a large oscilloscope read- 
ing error . The o profile of 
figure 12 illustrates this principle. 
The corrected a profile, using data 
recorded with a probe such that 
S = O.085, was computed on the basis of 
nine subregions whereas only four sub- 
regions were used for the profile of 
figure 10 (S = 0.23*0 because it was 
found that further subdivision intro- 
duced large uncertainty errors . It is 
also interesting, though not surpris- 
ing, that the smaller the value of 
y^/R, the smaller the resulting 
corrections . 

The other restriction concerns 
the choice of the radius, p n _ 1 of the 
innermost subregion. It is obvious 
that Pq_i must satisfy the inequality 
Pn-i > y^. Experience indicates that 
Qa-x should be chosen so that 
Qa-x > l.l4 y This restriction on 
the value of P n _ 1 places a lower limit 
of approximately 0.10 on the correc- 
tion factor ©£. Smaller values of 
this factor would magnify oscilloscope 
reading errors by an intolerable 
amount . 

It is difficult to estimate the 
overall accuracy of the preceding 
data-r eduction technique. The uncer- 
tainty resulting from the presence of 
random electromagnetic noise or 
oscilloscope reading errors has 
already been discussed. Another pos- 
sible error source is the neglect of 
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axial variations in a and U|| . However, the magnitude of such an error is 
probably small because the magnitudes of the theoretical current vectors, 

J^. and J^, decrease rapidly in the axial direction and because the effects of 

the higher upstream values of a and U|| may be canceled by the lower down- 
stream values . If the plasma stream is steady and fairly free from random 
electromagnetic noise, the data reduction technique presented in this section 
probably yields center-line values differing from the true values by no more 
than 20 percent . 

To test the validity of dividing the stream into several regions of con- 
stant cr, an experiment was performed using H 2 SO 4 solutions in plastic cylin- 
ders . The test was performed in the three concentric plexiglas cylinders as 


sketched in the plan view of figure 13 ; 


y 

! 




Plon view of plexiglas cylinders 



Figure 13.- Conductivity measured in acid 
solutions . 


the center-line position of the 
instrument is shown with the axis of 
the probe parallel to the z axis 
of the cylinders. The depth of each 
cylinder was at least l 8 cm because 
tests indicated that this value 
simulated a cylinder of infinite 
extent . All plexiglas walls were 
4.7 mm thick and the inside radii of 
the cylinders were 5-95 cm, 9*77 cm, 
and 12.3 cm as indicated in fig- 
ure 13 . The center cylinder was 
filled with a sulphuric acid solu- 
tion having a conductivity of 
74.8 mho/m as determined by a con- 
ductivity cell measurement. Simi- 
larly, the two annular spaces were 
filled with acid solutions having 
conductivities of 49 and 29 mho/m. 

It was found that the center- 
line value of could be pre- 

dicted from the theory given above 
if the probe were considered to be 
immersed simultaneously at the cen- 
ter of three cylinders of radii 
12.3 cm, 9-77 cm, and 5*95 cm con- 
taining electrolytes of conductivity 
29 j 20, and 25.8 mho/m, respectively. 
Then equation (47), 

— = 1585 mho/m-mV 

® E (S,P) 

could be used to compute the contri- 
bution of the currents induced in 
each cylinder to the total probe 
output;, for the instrument used in 
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this experiment, the conductivity calibration constant was 1585 mho/m-mV. 
Using the correction factors given in figure 5 and setting y^ = 2.22 cm, one 
may calculate the probe output as follows: 




1 

1585 



2.22 

, 12.3 



+ 20 © 2 ( 


2.22 
, 9 -77 



2.22 

5.95 



= TF5F [29 X 0.805 + 20 X 0.753 + 25-8 X 0 . 602 ] 
l>op 


= 0.0338 mV 


This theoretical value agrees quite well with the experimental value of 
0.034 mV. The off-axis experimental values of ^ for the case p / 0 could 
not be predicted by the theory because the two inner plexiglas walls created 
rather complicated boundary conditions . 


CONCLUDING REMARKS 


Tests in concentric cylinders filled with acid solutions furnished rea- 
sonable confirmation of the conductivity correction factors that were used to 
correct the raw conductivity data for the existence of a boundary and for 
cross-stream a variations. Using crU|| correction factors, the method was 
extended to raw crU|| data and the corrected a and aU^ profiles were used 
to compute a velocity profile. 

Although these results indicate that the present design has the advan- 
tages of practicality and theoretical justification, some additional develop- 
ment is needed. Specifically, the following items should be investigated. 

1. A new probe that permits variation of y^ should be constructed and 
tested for several values of 2p Q /y£ under identical plasma stream condi- 
tions and the corrected data from these tests should be analyzed to determine 
the significance of this parameter. It is believed that this parameter is a 
measure of the magnitude of mutual flow disturbance effects between adjacent 
support rods. 

^2. An oscillator having several output frequencies should be used with 
the probe to evaluate the importance of the frequency, to, of the impressed 
power . 


3. Improvement of the response time of the instrument system should be 
attempted . 

4 . Other materials, such as precision quartz tubing, should be tested 
for coil support rods. 
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5- A suitable experimental method for determining the oil., calibration 
constant Cy should be devised. 

6. The use of the voltages an(i $y P ro(iuce j means of elec- 
tronic division, a signal that is linearly proportional to U|j. An alternate 
probe system (primary and two fore and aft secondary coils located on the 
z axis; also, the axis of each coil alined with the z axis) should be 
investigated because a single support rod containing three coils would per- 
turb a plasma stream far less than a probe having three separate coil support 
rods . 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Nov. l4, 1967 
125-24-03-02-00-21 
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